Topological origin of the phase transition in a mean-field model 
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We argue that the phase transition in the mean-field XY model is related to a particular change 
in the topology of its configuration space. The nature of this topological change can be discussed on 
the basis of elementary Morse theory using the potential energy per particle V as a Morse function. 
The value of V where such a topological change occurs equals the thermodynamic value of V at 
the phase transition and the number of (Morse) critical points grows very fast with the number of 
particles N . Furthermore, as in statistical mechanics, also in topology the way the thermodynamic 
limit is taken is crucial. 
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It is customary in statistical mechanics to associate 
phase transitions with singularities in the equilibrium 
measure that describes the macroscopic system in its 
phase space. Recently it has been conjectured that 
second-order phase transitions might have a topologi- 
cal origin, i.e., that a thermodynamic transition might 
be related to a change in the topology of the configu- 
ration space, and that the observed singularities in the 
statistical-mechanical equilibrium measure and in the 
thermodynamic observables at the phase transition might 
be interpreted as a "shadow" of this major topological 
change that happens at a more basic level. Such a con- 
jecture has been put forward heuristically, based on nu- 
merical simulations where the averages and fluctuations 
— either time or statistical-mechanical — for some ob- 
servables of a geometric nature (e.g. configuration-space 
curvature fluctuations) related to the Riemannian ge- 
ometrization of the dynamics in configuration space have 
been computed. When plotted as a function of either the 
temperature or the energy, the fluctuations of the cur- 
vature have a singular behavior at the transition point 
which can be qualitatively reproduced using a geometric 
model. In such a model the origin of the singular behav- 
ior of the curvature fluctuations resides in a topological 
change [Q . Extensive numerical work || has shown that 
the curvature fluctuations indeed exhibit qualitatively 
the same singular behavior in many different models un- 
dergoing continuous phase transitions, namely (p 4 lattice 
models with discrete and continuous symmetries and XY 
models. The presence of a singularity in the statistical- 
mechanical averages as well as in the fluctuations of the 
curvature at the transition point has been proved analyt- 
ically for the mean-field XY- model || , which is the model 
that we will consider in the following. Moreover, a purely 
geometric, and thus still indirect, further indication that 
the topology of the configuration space might change at 
the phase transition has been obtained from numerical 



calculations for the ip 4 model on a two-dimensional lat- 

tice i-. 

In spite of all these indirect indications, no direct ev- 
idence, analytical nor numerical, of the presence of a 
change in the topology of the configuration manifold that 
can be related with a thermodynamic phase transition 
had yet been found. The purpose of the present Letter is 
to argue that such a topological change in configuration 
space indeed exists in the particular case of the mean- 
field XY model, and in doing so to clarify the nature of 
the relationship between topologic and thermodynamic 
transitions, discussing especially the role played by the 
thermodynamic limit. 

Let us now introduce the mathematical tools we need 
to detect and characterize topological changes, which are 
referred to as elementary critical point (Morse) theory 
|| . Such a theory links the topology of a given manifold 
with the properties of the critical points of functions de- 
fined on it. Given a (compact) manifold M and a smooth 
function / : M t— > R., a point x c S M is called a critical 
point of / if df — 0, while the value f(x c ) is called a crit- 
ical value. A level set / _1 (a) = {x e M : f(x) = a} of 
/ is called a critical level if a is a critical value of /, i.e., 
if there is at least one critical point x c <E / _1 (a). The 
function / is called a Morse function on M if its critical 
points are nondegenerate, i.e., if the Hessian of / at x c 
has only nonzero eigenvalues, so the critical points x c are 
isolated. Let us now consider a Hamiltonian dynamical 
system whose Hamiltonian is of the form 



1 N 

i=l 



V(<p) 



(1) 



where the <^'s and the 7r's are, respectively, the coordi- 
nates and the conjugate momenta. The dynamics of such 
a system is defined in the 2A r -dimensional phase space 
spanned by the </?'s and the 7r's, so that a natural choice 
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would be to investigate the topology of the phase space 
by using the Hamiltonian TC itself as a Morse function. 
However, for all the critical points of H., 7^ = Vi holds, 
so that there is no loss of information in considering as 
our manifold M the TV-dimensional configuration space, 
and the potential energy per particle V(<p) = V(<p)/N as 
our Morse function l| . Let us now consider the following 
family of submanifolds of M, 



M„ 



V -1 (-oo,«] = {ip S M : V(<p) < v} , 



(2) 



i.e., each M v is the set {^Pi}f =1 such that the potential 
energy per particle does not exceed a given value v. As 
v is increased from —00 to +00, this family covers suc- 
cessively the whole manifold M 0. All these manifolds 
have the same topology (homotopy type) until a critical 
level V _1 (w c ) = dM Vc is crossed; here the topology of M v 
changes. More precisely, the manifolds M v > and M v », 
with v' < v c < v", can not be smoothly mapped onto 
each other if v c is a critical value of V (if). A change in 
the topology of M v can only occur when v passes through 
a critical value of V. Thus in order to detect topological 
changes in M v we have to find the critical values of V, 
which means solving the equations 
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(3) 



Let us now apply this method to the case of the mean- 
field XY model. This model is particularly interesting 
four our purpose because the mean-field character of the 
interaction greatly simplifies the analysis || , allowing an 
analytical treatment of the Eqs. (^); moreover, a projec- 
tion of the configuration space onto a two-dimensional 
plane is possible. The mean-field XY model || describes 
a system of N equally coupled planar classical rotators. 
It is defined by a Hamiltonian of the class (Q) where the 
potential energy is 
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N 



2N 

Here ipi € [0, 2ir] is the rotation angle of the i-th rotator 
and h is an external field. Defining at each site i a clas- 
sical spin vector = (cos ipi, sinipi) the model describes 
a planar (XY) Heisenberg system with interactions of 
equal strength among all the spins. We consider only the 
ferromagnetic case J > 0; for the sake of simplicity, we 
set J = 1. The equilibrium statistical mechanics of this 
system is exactly described, in the thermodynamic limit, 
by mean-field theory |J . In the limit h —> + , the system 
has a continuous phase transition, with classical critical 
exponents, at T c = 1/2, or e c = 3/4, where e = E/N is 
the energy per particle. 

We aim at showing that this phase transition has its 
foundation in a basic topological change that occurs in 
the configuration space M of the system. Let us remark 



that since V(ip) is bounded, -h < V(tp) < 1/2 + h 2 /2, 
the manifold is empty as long as v < — h, and when v in- 
creases beyond 1/2 + h 2 /2 no changes in its topology can 
occur so that the manifold M v remains the same for any 
v > 1/2 + h 2 /2, and is an iV-torus. To detect topological 
changes we have to solve Eqs. (||). To this end it is use- 
ful to define the magnetization vector, i.e., the collective 
spin vector m = ■ir J2iLi s i: which as a function of the 
angles is given by 



/)= (^E cos ^^E sin ^) 

\ i=l i=\ / 



(5) 



Since, due to the mean-field character of the model, the 
potential energy (|4|) can be written as a function of m 
alone (remember that J = 1), the potential energy per 
particle reads 



V(tp) = V{m x ,m y ) = -(1- ml- ml) -hi 



(6) 



This allows us to write the Eqs. (||) in the form [i = 



1,...,A0 



(m x + h) sin tpi — m v cos ipi = 



(7) 



Now we can solve these equations and find all the critical 
values of V. The solutions of Eqs. (0) can be grouped in 
three classes: 

(i) The minimal energy configuration ifi = Vi, with a 
critical value v — Vq = —h, which tends to as h — > + . 
In this case, m 2 + m 2 = 1 . 

(ii) Configurations such that m y = 0, sin^j = 
These are the configurations in which ipi equals either 
or 7r; i.e., we have again ipi — Vz, but also the N con- 
figurations with ipk = 7r and tpi = V« 7^ k 7 then the 
N(N — 1) configurations with 2 angles equal to ir and 
all the others equal to 0, and so on, up to the configura- 
tion with ipi = 7r Vi. The critical values corresponding to 
these critical points depend only on the number of 7r's, 
n^, so that v(n„) = i[l - -^(A^-2^) 2 ] - j r {N~2n 7T ). 
We see that the largest critical value is, for N even, 
v{n v = N/2) = 1/2 and that the number of critical 
points corresponding to it is 0(2 N ). 

(iii) Configurations such that m x = — h and m y = 0, 
which correspond to the critical value v c = 1/2 + h 2 /2, 
which tends to 1/2 as h — > + . The number of these 
configurations grows with N not slower than N\ 0. 

Configurations (i) are the absolute minima of V, (iii) 
are the absolute maxima, and (ii) are all the other sta- 
tionary configurations of V. 

Since for v < v the manifold is empty, the topolog- 
ical change that occurs at Vq is the one corresponding 
to the "birth" of the manifold from the empty set; sub- 
sequently there are many topological changes at values 
vfon) € (vq,1/2] till at v c there is a final topological 
change which corresponds to the "completion" of the 
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manifold. We remark that the number of critical values 
in the interval [wo,l/2] grows with N and that eventu- 
ally the set of these critical values becomes dense in the 
limit N — > oo. However, the critical value v c remains 
isolated also in that limit. We observe that considering a 
nonzero external field h is necessary in order that V is a 
Morse function, because if h = all the critical points of 
classes (i) and (ii) are degenerate, in which case topologi- 
cal changes do not necessarily occur || . This degeneracy 
is due to the 0(2)-invariance of the potential energy in 
the absence of an external field. To be sure, for h 7^ 0, V 
may not be a Morse function on the whole of M either, 
but only on M v with v < v c , because the critical points 
of class (Hi) may also be degenerate, so that v c does not 
necessarily correspond to a topological change. 

However, this difficulty could be dealt with by using 
that the potential energy can be written in terms of the 
collective variables m x and m y — see Eq. (§). This im- 
plies that we consider the system of N spins projected 
onto the two-dimensional configuration space of the col- 
lective spin variables. According to the definition (|^) 
of m, the accessible configuration space is now not the 
whole plane, but only the disk 

D = {(m x ,m y ) :m 2 x +m 2 y < 1} . (8) 

Thus we want to study the topology of the submanifolds 

D v = {(m x ,m y ) G D : V{m x ,m y ) < v} . (9) 

The sequence of topological transformations undergone 
by D v can now be very simply determined in the limit 
h — > + (see Fig. |l]). As long as v < 0, D v is the empty 
set. The first topological change occurs at v = Vq = 
where the manifold appears as the circle m x + m y = 1 , 
i.e., the boundary 3D of the accessible region. Then as 
v grows D v is given by the conditions 

1 - 2v < m 2 x + m\ < 1 , (10) 

i.e., it is the ring with a hole centered in (0,0) (punc- 
tuated disk) comprised between the two circles of radii 
1 and V%v, respectively. As v continues to grow the 
hole shrinks and it is eventually completely filled as 
v = v c = 1/2, where the second topological change oc- 
curs. In this coarse-grained two-dimensional description 
in D, all the topological changes that occur in M between 
v = and v = 1/2 disappear, and only the two topologi- 
cal changes corresponding to the extrema of V, occurring 
at v — vq and v — v c , survive. This strongly suggests 
that the topological change at v c should be present also 
in the full TV-dimensional configuration space, so that 
the degeneracies mentioned above for the critical points 
of class (Hi) should not prevent a topological change. 

Now we want to argue that the topological change oc- 
curring at v c is related to the thermodynamic phase tran- 
sition of the mean-field XY model. Since the Hamiltonian 



is of the standard form ([!]), the temperature T, the en- 
ergy per particle e and the average potential energy per 
particle u = (V) obey, in the thermodynamic limit, the 
following equation: 

e = | + «(T), (11) 

where we have set Boltzmann's constant equal to 1. Sub- 
stituting the values of the critical energy per particle 
e c = 3/4 and of the critical temperature T c = 1/2 we 
get u c — u(T c ) — 1/2, so that the critical value of the po- 
tential energy per particle v c where the last topological 
change occurs equals the statistical-mechanical average 
value of the potential energy at the phase transition, 

v c = u c . (12) 

Thus although a topological change in M occurs at any 
N, and v c is independent of N, there is a connection of 
such a topological change and a thermodynamic phase 
transition only in the limit N —> 00, h —> + , when in- 
deed thermodynamic phase transitions can be defined. 
A similar kind of difference, as here between topological 
changes in mathematics (for all N) and phase transi- 
tions in physics (for N — » 00 only), also occurs in other 
contexts in statistical mechanics, e.g. in nonequilibrium 
stationary states [ p~T[ . 

The relevance of topologic concepts for phase transi- 
tion theory had already been proved, in the case of the 
two-dimensional Ising model Jl2| , though in a rather ab- 
stract context: in that case the phase transition was re- 
lated to a jump in the Atiyah index of a suitable vector 
bundle. However, Eq. (Q) strongly supports — albeit 
for a special model — the conjecture put forward in Ref. 
Q, i.e., that also the topology of the configuration space 
changes in correspondence with a thermodynamic phase 
transition. 

Since not all topological changes correspond to phase 
transitions, those that do correspond, remain to be de- 
termined to make the conjecture of Ref. jjj more precise. 
In this context, we consider one example where there are 
topological changes very similar to the ones of our model 
but no phase transitions, i.e., the one-dimensional XY 
model with nearest-neighbor interactions, whose Hamil- 
tonian is of the class ([j]) with interaction potential 

N N 

V(<p) = - ^ I 1 ~ cos^+i - <pi)] - h ^ cos fi ■ ( 13 ) 

i=l i=l 

In this case the configuration space M is still an TV-torus, 
and using again the specific interaction energy V — V/N 
as a Morse function we can prove that also here there are 
many topological changes in the submanifolds M v as v is 
varied in the interval [0, 1/2] (after taking h — > + ). The 
critical points are all of the type (fj = (pk = fl = • • • = TTj 
ipi = Vi 7^ j,k,l,...; however, at variance with the 
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mean-field XY model, it is no longer the number of 7r's 
that determines the value of V at the critical point, but 
rather the number of domain walls, rid, i.e., the num- 
ber of boundaries between "islands" of 7r's and "islands" 
of O's: v(nd) — na/2N. Since nd € [0,7V], the critical 
values lie in the same interval as in the case of the mean- 
field XY model. But now the maximum critical value 
v = 1/2, instead of corresponding to a huge number of 
critical points, which rapidly grows with N, corresponds 
to only two configurations with N domain walls, which 
are tp2k = 0, (fi2k+i = it, with k = 1, . . . , N/2, and the 
reversed one. 

Thus this example suggests the conjecture that a topo- 
logical change in the configuration space submanifolds 
M v occurring at a critical value v c is associated with a 
phase transition in the thermodynamic limit if the num- 
ber of critical points corresponding to the critical value v c 
is sufficiently rapidly growing with N. On the basis of the 
behavior of the mean-field XY model we expect that such 
a growth should be at least exponential. Furthermore, a 
relevant feature appears to be that v c remains an isolated 
critical value also in the limit N — > oo: in the mean-field 
XY model this holds only if the thermodynamic limit is 
taken before the h — > + limit: this appears as a topolog- 
ical counterpart of the non-commutativity of the limits 
h — > + and N — > oo in order to get a phase transition 
in statistical mechanics. 

We conclude with some comments. The sequence of 
topological changes occurring with growing V makes the 
configuration space larger and larger, till at v c the whole 
configuration space becomes fully accessible to the sys- 
tem through the last topological change. From a phys- 
ical point of view, this corresponds to the appearance 
of more and more disordered configurations as T grows, 
which ultimately lead to the phase transition at T c . We 
remark that the connection between the topology of the 
configuration space and the physics of continuous phase 
transitions made here via the potential energy, in partic- 
ular Eq. ([[4), only makes sense in the thermodynamic 
limit, where the potential energy per particle u(T) is 
well-defined since its fluctuations vanish then at least as 

Since a notion of universality arises quite naturally 
in a topological framework, it is tempting to think that 
universal quantities like critical exponents might have in 
general a topological counterpart. Finally, the fact that 
the topological changes appear at any iV opens a new 
possibility to study transitional phenomena in finite sys- 
tems, like atomic clusters, nuclei, polymers and proteins, 
or other biological systems. 
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FIG. 1. The sequence of topological changes undergone by 
the manifolds D v with increasing v in the limit h — ♦ + . 
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